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We report a density functional study of the electronic structure and magnetism of Y 2 Ni 7 . The 
results show itinerant magnetism very similar to that in the weak itinerant ferromagnet NiaAL The 
electropositive Y atoms in Y 2 Ni 7 serves to donate charge to the Ni host mostly in the form of s 
electrons. The non-spin-polarized state shows a high density of states at the Fermi level, N(Ef) 
due to flat bands. This leads to the ferromagnetic instability. However, there are also several much 
more dispersive bands crossing E{F), which should promote the conductivity. Spin fluctuation 
effects appear to be comparable to or weaker than NisAl, based on comparison with experimental 
data. Y 2 Ni 7 provides a uniaxial analogue to cubic NisAl for studying weak itinerant ferromagnetism, 
suggesting detailed measurements of its low temperature physical properties and spin fluctuations, 
as well experiments under pressure. 


I. INTRODUCTION 

Weak itinerant ferromagnetism is a topic of ongoing in¬ 
terest, both from the point of view of understanding the 
physical behavior of metals near quantum critical points, 
and because of the fact that these materials often have 
relatively high ordering temperatures when scaled to the 
ordered moment. There is also renewed interest in itin¬ 
erant magnetism because of the unusual m^netic prop¬ 
erties of the Fe-based superconductors, [ll-Q and spin- 
fluctuation pairing models for these and other unconven¬ 
tional superconductors. iHni 

Elemental fee Ni metal is a classic example of an itin¬ 
erant ferromagnet. Substitution by 25% with the triva- 
lent element A1 to form NisAl strongly reduces the Curie 
temperature to yield a weak itinerant ferromagnet near a 
critical point that can be reached under pressure. [illil 
Both NisAl and NisGa show evidence for strong ex¬ 
change enhancement and spin-fluctuation effects, includ¬ 
ing quantum spin fluctuation induced sup pre ssion of fer¬ 
romagnetism Bill in NisGa. [TgI - I^ These two 
compounds have very similar electronic structures, and 
differ mainly in the strength of the quantum spin fluctu¬ 
ations. This difference places them on opposite sides of 
a ferromagnetic quantum critical point at ambient pres¬ 
sure. Unusual physical behavior including non-Fermi liq¬ 
uid scalings extending to very low temperature has been 
found in NisAl under pressure. 0 

Y 2 Ni 7 forms in a rhombohedral (R3m) Gd 2 Co 7 struc¬ 
ture [23l-[^ and is a ferromagnet. [2614^ The mag¬ 
netism is unusual in that it has a high Gurie temperature 
of Tc ^ 54 K relative to the moment size, reported as 
m = 0.06 jiB - 0.08 jiB per Ni. [H, If one scales the 
Curie temperature of elemental Ni {Tc=627 K), which 
is also high relative to its moment size, by the square of 
the moment as usual, one would infer an expected Curie 
temperature of only ^ 6-7 K for Y 2 Ni 7 . The system is 
also unusual in that in spite of the high Tc, it is very sen¬ 
sitive to alloying, both by H incorporation, and by 
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metal alloying. IHIsl Here we investigate the magnetic 
and electronic properties of Y 2 Ni 7 in relation to Ni metal 
as well as the weak itinerant ferromagnet NisAl and the 
related compound NisGa. We find behavior reminiscent 
of NisAl. 


II. METHODS AND STRUCTURE 

The present density functional calculations were done 
using the generalized gradient approximation of Per dew, 
Burke and Ernzerhof and the linearized augmented 
planewave (LAPW) method as implemented in the 
WIEN2k code. We also performed calculations us¬ 
ing the local spin density approximation (LSDA) and we 
also tested the effects of spin-orbit coupling. We used 
well converged basis sets, with a planewave cutoff, Kmax 
determined by RminKmax=^i where Rmin is the small¬ 
est sphere radius, here 2.25 Bohr for both Y and Ni. The 
calculations were done using the experimental lattice pa¬ 
rameters, [ill, a=4.947 A, c=36.25 A. The internal co¬ 
ordinates were determined by total energy minimization. 
Eor this purpose we started with the structure of the pro¬ 
totype (Gd 2 Co 7 ). The resulting structure is given in Ta¬ 
ble Hand depicted in Eig. [H Plasma frequencies were de¬ 
termined using the optical package of the WIEN2k code, 
which uses integration of the squared band velocity on 
the Eermi surface for this purpose. In this code the ve¬ 
locities come from calculations of the dipole (momentum) 
operator. 


III. RESULTS AND DISCUSSION 

Y 2 Ni 7 , which is ^78% Ni, has metal atoms in dis¬ 
torted 12-fold cages. As mentioned, NisAl, which con¬ 
tains 75% Ni, also with a trivalent element has suppressed 
ferromagnetism relative to Ni and is near an interesting 
quantum critical point. [H, [l9|, [22| NisAl has Tc=41.5 
K, M=0.08 /i^/Ni, similar to Y 2 Ni 7 . One signature of 
physical importance of fluctuations associated with the 
quantum critical point in NisAl is an overestimate of the 
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FIG. 3. Spin density plot showing the orbital character of the 
magnetization. The left panel shows a view along the c-axis 
while the right panel is perpendicular to the c-axis. The labels 
indicate the Ni atoms each layer along c, as in Fig. [T] Note 
the pronounced anisotropy around the Ni atoms. 


FIG. 1. Grystal structure of Y 2 Ni 7 showing the layering of 
the metal atoms. Y is shown as light grey, with Ni as dark 
blue. The numbers denote the atom number as in Table [H 



FIG. 2. Band structure of non-spin-polarized Y 2 Ni 7 as ob¬ 
tained with the PBE GGA in a scalar relativistic approxima¬ 
tion (left) and with spin-orbit coupling (right). 


ordered moment in standard density functional calcula¬ 
tions. In the case of NisAl, the calculated spin moment 
in the local density approximation is /i^/Ni. 

M 

For Y 2 Ni 7 we obtain a spin magnetization of 1.29 
per formula unit (f.u.) with the PBE GGA, with a mag¬ 
netic energy of 0.04 eV/f.u. On a per Ni basis, this is 
0.18 /iB and 5.7 meV ^ 70 K, i.e. only slightly higher 
than Tc'=54 K. This is indicative of being in the strongly 
itinerant (Stoner) limit. It is notable that the moment 
size is ^ 0.1 fiB higher than the reported experimental 
value. The spin density is illustrated in Fig. [3l The mo¬ 
ments on the different Ni sites, as determined by the spin 
magnetization in the corresponding LAPW spheres, vary 
considerably as seen in Table [Tl There is a small back po¬ 
larization in the interstitial and around the Y atoms, as 
is commonly found in 3d transition metal ferromagnets. 
This interstitial spin density is derived from extended or¬ 
bitals, such as the metal s states. The two Ni sites with 
low numbers of Ni neighbors (Nil and Ni5) have the low¬ 


est moments, with Ni5 by far the lowest. However, there 
is not a clear correlation between coordination and the 
Ni moment for the other sites. 

Calculations for fee Ni done in the same way yield an 
ordered moment of 0.635 /ib and a magnetic energy of 
0.062 eV/atom ^740 K (c.f. Tc'=627 K). Spin orbit cou¬ 
pling has only a very small effect, reducing the spin mo¬ 
ment by less than 0.01 /iB/hu. and producing orbital 
moments of 0.007 iib - 0.025 i^b per Ni, depending on 
the specific site, always aligned with the spin moment, 
as expected from the third Hand’s rule. The effect of 
spin orbit coupling on the electronic structure is also very 
small. This is seen in the band structure near the Fermi 
energy, which is given in Fig. [2] in a scalar relativis¬ 
tic approximation and with spin orbit coupling. In the 
following, we give scalar relativistic results. 

With the LSDA we still obtain an overestimation of 
the ordered moment, although somewhat smaller, specif¬ 
ically, Mlsda=^-^7 /iB/hu. and SElsda=0.027 eV/f.u. 
or ^47 K on a per Ni basis, i.e. slightly less than T^. This 
similarity of the magnetic energy to the experimental Tq 
is similar to what was found for NisAl, suggesting 
a similarity of Y 2 Ni 7 and NisAl. However, it should be 
noted that the more complicated non-cubic structure of 
Y 2 Ni 7 may lead to larger experimental uncertainty in the 
determination of the moment size due to magnetocrys¬ 
talline anisotropy and more possibilities for intrinsic de¬ 
fects. 

The calculated electronic density of states and projec¬ 
tions of Ni d and Y d character are given in Fig. IH both 
for the non-spin-polarized and the ferromagnetic cases. 
As shown, there is a sharp peak in the density of states 
almost exactly dX Ep- The high value leads to a Stoner 
instability and ferromagnetism, with an exchange split¬ 
ting of the Ni d bands, although not in a perfect rigid 
band fashion (note the change in the shape of the ex¬ 
change split peak between majority and minority spin). 
Also, the Y d states are above Ep in this compound. In 
contrast, Y metal is a transition element with ^ 2 d elec- 
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TABLE L Calculated internal structural parameters and 
atomic moments with the PBE GGA of Y 2 Ni 7 at the exper¬ 
imental lattice parameters of a=4.947 A, c=36.25 A, space 
group 166, R3m with hexagonal coordinates. The nearest 12 
neighbors are given as “coord.”. 


Atom 

X 

y 

z 

coord. 

m(/is) 

Y1 6c 

0 

0 

0.0504 

12 Ni 

-0.03 

Y2 6c 

0 

0 

0.1472 

12 Ni 

-0.02 

Nil 18h 0.1671 

0.8329 

0.4430 

7 Ni, 5 Y 

0.17 

Ni2 9e 

1/6 

1/3 

1/3 

8 Ni, 4 Y 

0.33 

Ni3 6c 

0 

0 

0.2782 

9 Ni, 3 Y 

0.22 

Ni4 6c 

0 

0 

0.3883 

9 Ni, 3 Y 

0.23 

Ni5 3b 

0 

0 

1/2 

6 Ni, 6 Y 

0.10 


trons. This implies a charge transfer from the Y d states 
to the Ni host matrix. This is not surprising in view 
of the electropositive nature of Y relative to Ni. Based 
on integration of the Ni d density of states in compari¬ 
son with fee Ni, most of this charge transfer is to the s 
electrons. This inference is based on the fact that the d 
electron count of Ni does not show an increase commen¬ 
surate with the additional charge. 

Within Stoner theory [36|, [13 the susceptibility of a 
metal is given by a random phase approximation (RPA) 
formula, y = xo/[^ — N{Ef)I]^ where xo is the bare Pauli 
susceptibility, xo = with appropriate units. 

This formula is exact at the level of band structure calcu¬ 
lations, but neglects the effects of spin-fluctuations which 
can renormalize the spin susceptibility (see Ref. [HI for a 
detailed discussion applied to Pd, which is a high suscep¬ 
tibility paramagnetic metal). The Stoner theory itiner¬ 
ant ferromagnetic instability occurs when N{Ef) = ^ 

which is the point where the RPA enhancement factor 
1/[1 — N{Ef)I] diverges. 

In contrast to standard local moment magnetic mate¬ 
rials, in the itinerant limit a metal with properties de¬ 
termined by the non-spin-polarized electronic structure 
occurs above Tc- There is typically a second order or 
near second order phase transition at Tq. In Y 2 Ni 7 this 
is a rather interesting metallic state in part because of 
the high density of states. 

We obtain N{Ef)=‘24.0 eV“^ per f.u. on a both spins 
basis. This well above the criterion for Stoner mag¬ 
netism, B S [13 ^ where I is 0.7-0.9 eV for 

heavy 3d transition elements [HI and N is N{Ef) ex¬ 
pressed on a per atom per spin basis (1.71 eV“^ in the 
present case). The fixed spin moment energy (Fig. [5)) as 
a function of magnetization therefore drops quickly from 
zero but then rises because of the narrowness of the peak. 
The curve is featureless and smooth except for the min¬ 
imum corresponding to the ferromagnetic ground state. 
There is no sign of any metamagnetic state at higher 
magnetization. 

The peak at Ef is derived from Ni d orbitals on the 
Ni2, Ni3 and Ni4 sites, with a somewhat smaller contri¬ 
bution from the Nil site and practically no contribution 




FIG. 4. Electronic density of states of Y 2 Ni 7 as obtained with 
the PBE GGA for non-spin-polarized (top) and ferromagnetic 
(bottom) states. 



m (lae/f-u ) 


FIG. 5. Fixed spin moment energy as a function of con¬ 
strained magnetization. The inset shows an expanded scale 
for low magnetizations. 
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from the Ni5 site. This mirrors the distribution of the 
moments in the ferromagnetic state. For an atom with 
trigonal site symmetry (i.e. axial with a three fold axis), 
as are the Ni3, Ni4 and Ni5 atoms, the d electron crystal 
field levels are a singly degenerate ag with ^ along the 
axis), and two double degenerate eg levels — y^^xy 
and xz-\-yz). The Ni3 and Ni4 contributions to the peak 
at Ep are from the Cg {xz-\-yz) orbital. The Nil and Ni2 
contributions are from the same two d orbitals in this ref¬ 
erence frame, but not in equal proportion reflecting the 
lower site symmetry. The orbital character is reflected in 
the spin density (Fig. [3]), where the lobes corresponding 
to these orbitals are seen around the Ni sites. As seen 
in the figure these orbitals are not oriented favorable for 
bonding interactions, which is a fact consistent with the 
narrowness of the peak. 

The magnetic behavior found in the calculations is con¬ 
sistent with what is expected from extended Stoner the¬ 
ory based on the sharply peaked density of states. Within 
extended Stoner theory, [13, |4l| which involves a rigid 
band approximation, one exchange splits the density of 
states to obtain ma gnetiza tion. The stat ionary solutions 
are points for which N{m) = where N{m) is the av¬ 
erage density of states between the position of the Fermi 
level for minority spin and that for majority spin to ob¬ 
tain the magnetization, m using the non-spin-polarized 
density of states. 

In other words N(m) is the magnetization divided by 
the energy shift between majority and minority spin den¬ 
sities of states needed to produce this magnetization. 
Thus a sharp peak with little area underneath it will 
yield a strong initial ferromagnetic instability in a fixed 
spin moment plot, but will not lead to a high magnetiza¬ 
tion. This is the origin of the high Tc with low moment. 
In other words a narrow peak with a very high N(Ef) 
leads to a strong initial instability, i.e. a large negative 
X = Xo/[l ~ E{Ef)I]- However if the weight under the 
peak i s small, the net magnetization will be low because 
N{m) will not be large for finite m since it is an integral. 
Note also that the extended Stoner formula uses the con¬ 
nection between the moment size and the exchange split¬ 
ting through the fact that the moment is an integral of 
the density of states over the energy range coming from 
the exchange splitting. This connects the moment, the 
exchange splitting and the magnetic energy. 

The density of states at the Fermi level is reduced to 
N{Ef)=9.1 eV“^ per f.u., mainly from the minority spin 
in the ferromagnetic ground state. The majority and 
minority spin contributions are N^{Ef)= 2.8 eV“^ and 
Ni{Ef)=^-S eV“^, respectively. This corresponds to a 
bare specific heat coefficient, 76 are= 18.9 mJ/(mol K^) on 
a per f.u. basis. It would be interesting to compare with 
experiment to determine the specific heat enhancement, 
717bare^ which if large might be an indicator of quantum 
spin fluctuations. 

We now return to the non-spin-polarized case, which 
should correspond to the electronic structure above Tc 
in this itinerant material. The band structure (Fig. [2]) 


shows heavy bands close to Ep, specifically the band 
just below Ep along a large part of the T-L line and 
the band at Ep along P-Z. In addition there are sev¬ 
eral much more dispersive bands crossing Ep- The con¬ 
sequence is that although the value of N{Ef) is high, 
the material can have a reasonable conductivity. This 
multi-sheet Fermi surface characteristic also occurs in 
NisGa and NigAl. [H HH In general the conductiv¬ 
ity can be written as a oc cj^r, where r is an effec¬ 
tive inverse scattering rate and uOp is the plasma fre¬ 
quency. We obtain plasma energies kujp of Qp^a=‘^-^ eV 
and flp^c=2.0 eV, for the basal plane and c-axis direc¬ 
tions, respectively. The implied conductivity anisotropy 
is modest, cFalcfc ^1-3. The anisotropy increases in the 
low temperature ferromagnetic state, for which we ob¬ 
tain Gajcfc ^ 1 . 8 . Also, even though the magnetization is 
small we obtain a significant transport spin polarization, 
defined by P = (cr^ — cr^)/(cr| + cr^). We obtain P=0.09 
both in the basal plane and c-axis directions. 

As mentioned, Y 2 N 7 has a remarkably high ordering 
temperature in relation to its moment. Most magnetic 
materials are described in terms of local moments and 
their interactions through the interatomic exhange cou¬ 
plings, These describe spin wave dispersions, which 
are transverse in character. Longitudinal degrees of free¬ 
dom, which correspond to changes in the local moment 
size are hard and not involved in the phase transition. 
This leads to a simplification in which one can treat 
the phase transition using the effective spin Hamilito- 
nian and the effect of increasing temperature as the ex¬ 
citation of spin-waves. This local moment case has been 
well described theoretically and numerical simulations 
of the temperature dependent magnetic properties and 
phuse transitions are practical even for complex systems. 
1^ . 1^ Importantly, the paramagnetic state above the 
ordering temperature is a disordered local moment state 
in which the local moment directions may be regarded as 
fluctuating in time but retaining their size. This means 
that the atomic Hund’s rule energy associated with the 
moment formation (i.e. the longitudinal degree of free¬ 
dom) is not involved in the phase transition as this con¬ 
tribution to the magnetic energy is present in both the 
ordered and paramagnetic phases. 

The itinerant limit has excitation of both transverse 
and longitudinal degrees of freedom with temperature 
through coupling to the electronic system and is not de- 
scribable by an effective spin Hamiltonian. This is more 
difficult to treat theoretically as it involves coupling to 
the electrons without separation of electronic and mag¬ 
netic degrees of freedom. ¥ 2^7 and NisAl are both ap¬ 
parently close to this limit, and based on comparison of 
density functional results with experiment both also have 
renormalizations of their ground states due to quantum 
fluctuations. Itinerant magnets also have magnetic con¬ 
tributions to the energy in the paramagnetic state, Q 
above but close to the ordering temperature. However, 
these are reduced as the moment size is reduced, and in 
the itinerant limit become negligible. This means that 
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in the itinerant limit all of the magnetic energy includ¬ 
ing the onsite Hand’s energy is available to drive the 
ordering, providing an explanation for the high ordering 
temperatures relative to the moments in these itinerant 
materials. In other words, in itinerant magnets disorder¬ 
ing implies destroying the moments, which has an energy 
cost from the on-site Hand’s energy and this leads to high 
ordering temperatures. We note that the Stoner model 
of itinerant magnetism has a parallel in the Slater model 
of itinerant antiferromagnetism and that high ordering 
temperatures in certain antiferromagnets have been dis¬ 
cussed in a way similar to the above. In neither 

case (itinerant or local moment magnets) can energy dif¬ 
ferences by themselves be simply interpreted as the or¬ 
dering temperature. 


IV. SUMMARY AND CONCLUSIONS 

We report density functional calculations of the elec¬ 
tronic structure and magnetic properties of ¥ 2 ^ 7 . 
Y 2 Ni 7 shows similarity to the weak itinerant ferromagnet 
NisAl, which has a modest renormalization of the mag¬ 


netism due to nearness to a quantum critical point. The 
overestimation of the magnetization relative to experi¬ 
ment is a bit smaller in Y 2 Ni 7 implying perhaps some¬ 
what weaker fluctuation effects. This is in contrast to 
NisGa, where magnetic ordering is apparently completely 
suppressed by quantum spin fluctuations. The weak itin¬ 
erant ferromagnetism in Y 2 Ni 7 arises from a Stoner insta¬ 
bility of a rather interesting metallic state. This state fea¬ 
tures a very narrow density of states peak at a mix¬ 
ture of dispersive and flat bands crossing Ef and modest 
but non-negligible anisotropy of the plasma frequency. 
Thus Y 2 Ni 7 is a uniaxial analogue of cubic NisAl. It will 
be of interest to study the low temperature properties 
of Y 2 Ni 7 , its magnetic fluctuations and the pressure de¬ 
pendence of the magnetic, thermodynamic and transport 
properties in comparison with NisAl. 
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